home *** CD-ROM | disk | FTP | other *** search
/ Power Programmierung 2 / Power-Programmierung CD 2 (Tewi)(1994).iso / c / library / dos / diverses / leda / src / plane / _point.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-11-15  |  2.7 KB  |  123 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  2.1.1                                                 11-15-1991
  4. +
  5. +
  6. +  _point.c
  7. +
  8. +
  9. +  Copyright (c) 1991  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #include <LEDA/plane.h>
  17. #include <math.h>
  18. #include <ctype.h>
  19.  
  20. static const double eps = 1e-10;
  21.  
  22.  
  23. //------------------------------------------------------------------------------
  24. // points 
  25. //------------------------------------------------------------------------------
  26.  
  27. point_rep::point_rep()  { count=1; x = y = 0.0; }
  28.  
  29. point_rep::point_rep(double a, double b) 
  30. { x = a; 
  31.   y = b; 
  32.   count = 1; 
  33. }
  34.  
  35.  
  36.  
  37. point::point()                  { ptr = new point_rep; }
  38. point::point(int)               { ptr = new point_rep; }
  39. point::point(ent x)             { ptr= (point_rep*)x; ptr->count++;  }
  40. point::point(double x, double y){ ptr = new point_rep(x,y); }
  41. point::point(vector v)          { ptr = new point_rep(v[0], v[1]); }
  42. point::point(point& p)          { ptr = p.ptr; ptr->count++; }
  43. void point::clear()             { if (--(ptr->count)==0) delete ptr; }
  44.  
  45. point point::rotate(point origin, double alpha)
  46. { if (origin == *this) return *this;
  47.   segment s(origin,*this);
  48.   return s.rotate(alpha).end();
  49. }
  50.  
  51. point point::rotate(double alpha)
  52. { return rotate(point(0,0),alpha);
  53.  }
  54.  
  55. point point::translate(double alpha, double d)
  56. { double dx = cos(alpha) * d;
  57.   double dy = sin(alpha) * d;
  58.   return point(ptr->x+dx,ptr->y+dy);
  59.  }
  60.  
  61. point point::translate(const vector& v)
  62. { return point(ptr->x+v[0],ptr->y+v[1]);
  63.  }
  64.  
  65. double point::distance(point p) 
  66. { return hypot(p.ptr->x - ptr->x, p.ptr->y - ptr->y); } 
  67.  
  68. double point::distance()        
  69. { return distance(point(0,0)); }
  70.  
  71. point& point::operator=(const point& p)
  72. { p.ptr->count++;
  73.   if (--ptr->count == 0)  delete ptr;
  74.   ptr = p.ptr;
  75.   return *this;
  76. }
  77.  
  78. int point::operator==(const point& p) 
  79. { return (fabs(ptr->x - p.ptr->x) < eps && fabs(ptr->y - p.ptr->y) < eps); }
  80.  
  81.  
  82. ostream& operator<<(ostream& out, const point& p)
  83. { out << "(" << p.xcoord() << "," << p.ycoord() << ")"; 
  84.   return out;
  85.  } 
  86.  
  87. istream& operator>>(istream& in, point& p) 
  88. { // syntax: {(} x {,} y {)}
  89.  
  90.   double x,y; 
  91.   char c;
  92.  
  93.   do in.get(c); while (isspace(c));
  94.   if (c != '(') in.putback(c);
  95.  
  96.   in >> x;
  97.  
  98.   do in.get(c); while (isspace(c));
  99.   if (c != ',') in.putback(c);
  100.  
  101.   in >> y; 
  102.  
  103.   do in.get(c); while (c == ' ');
  104.   if (c != ')') in.putback(c);
  105.  
  106.   p = point(x,y); 
  107.   return in; 
  108.  
  109.  } 
  110.  
  111. int compare(point& a, point& b)
  112. { 
  113.   real xa = a.ptr->x;
  114.   real xb = b.ptr->x;
  115.   real ya = a.ptr->y;
  116.   real yb = b.ptr->y;
  117.   int r = compare(xa,xb);
  118.   if (r==0) return compare(ya,yb);
  119.   return r;
  120. }
  121.  
  122.  
  123.